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We study the short-time dynamics of a liquid ligament, held between two solid cylinders, when one 
is impulsively accelerated along its axis. A set of one-dimensional equations in the slender-slope 
approximation is used to describe the dynamics, including surface tension and viscous effects. An 
exact self-similar solution to the linearized equations is successfully compared to experiments made 
with millimetric ligaments. Another non-linear self-similar solution of the full set of equations is 
found numerically. Both the linear and non-linear solutions show that the axial depth at which the 
liquid is affected by the motion of the cylinder scales like \/t. The non-linear solution presents the 
peculiar feature that there exists a maximum driving velocity U'^ above which the solution disap¬ 
pears, a phenomenon probably related to the de-pinning of the contact line observed in experiments 
for large pulling velocities. 


When the surface of a cohesive material is moved impulsively, for instance because it is impacted or 
pulled by another, the core of the material develops internal stresses and deformations. If the material is 
incompressible like liquids when moderately shocked, it is its boundary itself which must deform. For a 
solicitation localized in space or time (or both), raises the natural question of the extent and dynamics of 
the surface perturbation. Classical examples with solids include the celebrated Hertz contact problem, and 
the Saint-Venant solution for a bar impacted longitudinally, of which Hopkinson gave a variant when it is 
pulled (see Lov^and GraflPfor historical references and a detailed exposition). In this precise problem, the 
stress is transported along the bar at the velocity of sound, according to a non-dispersive wave equation so 
that the deformation at any location in the bar reflects that at the impact coordinate an instant earlier. 

Waves, however, propagate at the surface of liquids at different speeds according to their wavelength. 
Kelvin and RayleiglP, studying the pattern of standing waves at the surface of running water past a fixed 
obstacle have shown that their wavelength is short, capillary dominated upstream of the obstacle, and that 
longer wavelengths, gravity dominated, are found downstream of it: dispersive waves are sorted depending 
on the direction they propagate. 

While the stress penetration into a solid bar allows solving for its stability, and ultimate breakupl^, the 
dispersive dynamics of capillary waves on a liquid ligament helps understanding the remnant liquid mass 
attached to a rod quickly removed from a pooP. Impacts distort l ic j uid surfaces into a variety of non-trivial 
shape^^M^, in particular when they are mediated by capillary wave^®^ Here, we provide an original solution 
for the short time response of a liquid ligament of which one extremity is either pulled, or pushed, and in 
all cases suddenly accelerated axially. 

We consider a liquid ligament, initially cylindrical of length and radius R, held between two cylinders 
of same radius. The initial state is then characterized by the aspect ratio 2o/R and the Ohnesorge number 
Oh = r]/y/paR where p is the density of the fluid, a its surface tension and p its dynamic viscosity. One of the 
cylinder is impulsively accelerated at t = 0 along its axis with a velocity V, either positive or negative (see 
figure Ih. The strength of this impulsion is measured by the dimensionless parameter Uq = V/y^a/pR. In 
the foll Qwing, length and time scales are made dimensionless using the radius R, and the capillary timescale 
Tc = ^/pWJa. 

We describe the dynamics of the ligament by a one-dimensional model 

{h% = -{hy,, ( 1 ) 

Ut = -uu,+K, + 30h3^^, ( 2 ) 

where 

K = __ ^ 

(1 + /i2)3/2 

is the full curvature of the axisymmetric interface, u{z, t) is the local axial velocity, h{z, t) is the local radius, 
z is the axial coordinate, and t is the time variable. Although these equations have been derived within the 
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FIG. 1. Short-time dynamics of the stretched ligament for Oh=10“^ and Uq — —2. The red dashed curves correspond 
to the velocity in the ligament. The graduations are the same in both directions and the unit distance corresponds 
to one cylinder radius. A(t) is the characteristic depth at which the velocity has decreased significantly and is clearly 
much larger than the cylinder displacement. 


slender-slope approximatioiP^^, the y h ave proven to be remarkably efficient at describing the jet dynamics 
when the slope is far from being smalP^, and even close to the ligament breakup, where the slope diverge^^^. 
Also, choosing the full curvature (§ instead of the asymptotical limit for small slopes guarantees that the 
static equilibrium shapes can be correctly recovered. This one-dimensional idealization is expected to be 
relevant to describe the flow in the vicinity of the moving wall at the contact of the rod in 2 ; = 0 where the 
velocity is, by cohesion, enforced to be uniform, equal to the rod velocity. Indeed, we assume in the present 
analysis that the contact line is fixed to the boundary, an hypothesis justified for small values of [/q, but 
which fails for larger values, as we will see. 

Linearizing equations 0 and 0 for small velocities Uq and small deformations around ho = 1 gives for 
the velocity field 


2uu + Uzz + Uzzzz - 6 Ohuzzt = 0, (4) 

from which one deduces the dispersion relatiorP^ 

— 2ci;^ — ^ — 6i Oh k^u = 0, (5) 

which connects the frequency uj to the wavenumber k of a normal mode The radius h is related to 

u through 2ht = —Uz. As we expect the dynamics in the first instants to be localized close to the moving 
boundary, we assume that it is not influenced by the finite size £o/R of the ligament. More precisely, we 
consider that the other boundary is sufficiently far so that we can make lo/Ft 00 in the analysis, an 
hypothesis on which we come back below. 

It is convenient to solve the problem in the frame moving with the boundary which experiences an 
acceleration UoS{t), i.e. a velocity jump Uo'H{t), where 6{t) and 'H(t) are the Dirac delta function and the 
Heaviside step function respectively. The acceleration a = dfU in this non-galilean reference frame can then 
be written as 


a{z^t) = aA{z,t) — UoS{t) =—^ / a{uj,z)e — UoS{t)^ ( 6 ) 

J_oo 

where the boundary conditions on the absolute acceleration aA{z^t) give for d 

d(ci;,0) = l; d(co’,oc) = 0; dza{uj,0) = 0; dza{uj, oc) = 0. (7) 

Note that the conditions on dzh comes from the fixed contact line condition which requires dzU to vanish 
on the boundary. Since a a satisfies we obtain d satisfying Q as 

ik+(uj)z 


a{uj, z) = C-^{uj)e' 


( 8 ) 
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c+ = - 




where k±{uj) are the two solutions of ^ such that ^m{k±) > 0 


k±{uj) 



Qi Oh u) ± ■\/(l + 6i Oha;)2 + 


(9) 


( 10 ) 


The functions k± can be analytically defined on the integration contour if it is a line parallel to the ujr 
axis above the maximum growth rate, that is above, say, uji = 1/2. This choice of contour also guarantees 
causality: for negative times, the integration contour is closed in the upper plane such that the solution is 
null for t < 0. From the expression of a, we deduce the expressions of u and h 


u{z, t) 


Uo f 

27r 


—iuj 




(11) 


h{z,t) = 1 + 


Uo f 

47 r J_ 






( 12 ) 


For small times, we expect the solution to vary close to the origin only. In this region, the main contribution 
to the integrals is then expected to come from large values of |cc;| for which \k±{uj)\ is also large and given 
by 


k± ^ Aj 


with 


A| = 3 Oh±y9 Oh^-2. 


(13) 


The functions A± are both real and positive for Oh > \/2/3, or complex conjugate with real(A±) > 0 for 
Oh < \/2/3. The coefficients c+ and c_ become independent of uj. The above expression for a, u and h can 
then be evaluated for t > 0 using known contour integrals, and we obtain 


a{z,t) = 


Uoz 


2^/7^t^/i - 


C+A+ + c_A_ 


(14) 


i{z,t) = -Uo 




c_Ert(^ 

\2Vi 


(15) 


/i(^, t) — 1 “h U^^/t 


c^- 




At 


c+\% 


2Vi 


-Evk (^ 

\2Vi 


c_A_ - c-X_z^ ^ 

-4t-^Erfc i ^ , 

2^/t \2^/t J 




(16) 


Note that the above expressions are the sum of complex conjugate terms when Oh < \/2/3, while all the 
terms are real when Oh > \/2/3. As soon as Oh is nonzero, u and h tend exponentially fast to a constant 
on a length scale proportional to ^JTfOh. For small time, the dynamics is therefore limited to the close 
neighborhood of the moving boundary, as assumed. More precisely, the finite size of the ligament is not 
expected to affect the dynamics if t remains smaller than Oh(^o/^)^- 

It is worth mentioning that expression (16) for h satisfies volume conservation: The ligament volume 
remains constant up to 0{Uo) terms since SV = tt J^{h‘^{z) — l)dz ^ TrUot balances the volume —irUot 
associated with the elongation {Uq < 0) or compression {Uq > 0) of the ligament. Expression (16) for h can 
be written in the self-similar form 


h{zU) = 1 + UoVi X H{z/Vt). 


( 17 ) 
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FIG. 2. Linear and non-linear self-similar solutions for Oh = 0.0025. Since Uq appears in equations (21) and (22), 
elongated and compressed shapes are symmetric with respect to the unperturbed interface. 


The scaling variable z/^/i reflects the form of the dispersion relation ^ which is tangent to ^ at 
sho rt time and small scale, a feature shared, in particular, with the dynamics of flexural waves on solid 
rod^SHmSl the result in ([T^ is besides reminiscent of those established by Boussinesq for vibrating elastic 
beamJ^. It is genuine to the imposed transverse length scale i?, and differs from the one known to prevail 
in surface-tension-driven flows ^ in absence of externally imposed length scal^^^. Yet, the scaling law 
in ^/i can be obtained using the simple ansatz derived from the boundary conditions 


u{z, i) = C/o f/ (I) , h{z, t) = 1 + (I) , 


(18) 


where A(t) and /i(t) are the length scale defined in figure This ansatz assumes that velocity and surface 
distorsions follow the same scaling, an hypothesis which is commoiP, but not systematicP. The amplitude 
scaling for the velocity directly comes from the boundary condition which prescribes that u{z^ t) should be 
equal to —Uq at infinity for all time. Plugging (18) in equations 0 and ([^ and requiring the balance of all 
linear terms, we get 


/i(t)A'(t) ^ /i'(t)A(t) ^ Uq and /i(t) ^ X{t)Uo, 


(19) 


giving 


= Vi and /i(t) = UoVi. 


( 20 ) 


Expression (16) for h has been obtained in the limit t ^ 0 and \Uo\ 1. A more general self-similar 
solution of equations 0 and 0 can be obtained by removing the constraint |I7o| ^ 1- Indeed, inserting the 
ansatz (18) with ( [2Q| ) in equations 0 and we obtain, for small times, a system of non-linear equations 
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FIG. 3. Slope of the interface at the contact line {—hz{0,t) = cotan(^)) as a function of the driving velocity Uq, for 
Oh = 10~^. The limiting velocity is U'^ ~ 2.2. The pictures in inset correspond to experiments for which the driving 
velocity is smaller than f/* (left: Uq ~ 1.17) and bigger than f/* (right: Uq ~ 3.73) : in the later case, the contact 
line is not pinned anymore and dewetting occurs. The red dashed angle drawn on the first figure corresponds to the 
theoretical slope. 


with respect to the self-similar variable x = z/\/t ^OT the function U and H 


H{x)-xH\x) = -U\x), 


--xU'{x) 


H"'{x) 


3U^H"^{x)H'{x) 


(1 + U^H'^ix))^/^ (1 + U^H'^ix))^/^ 


3 0hU"{x). 


( 21 ) 

( 22 ) 


Note that the small t hypothesis guarantees that nonlinearity only appears in the curvature term. This 
system of ordinary differential equations can be solved numerically. The result of the integration is displayed 
in figure for different values of Uq and Oh = 0.0025. The nonlinear solution is very close to the linear 
solution as soon as T/q < 1, a property that we have also observed for other values of Oh. 

An interesting feature of equations (21) and (22) arises when seeking a solution for large values of the 


driving velocity Uq. Indeed, there exists a limit value of Uq above which there is no solution to this set of 
equations. Figure shows the slope of the interface at the contact line, as a function of the driving velocity 
Uq for Oh = 0.01. For Uq < the linear behavior described by equation ( p!^ is in agreement with the 
numerical results : the slope of the interface is well approximated by l/o/\/2, which is the value of —hz{0, t) 
given by equation (16) for Oh = 0. 

Above T/q = 1, the solution departs from the linear solution, and finally diverges for a finite value of the 
velocity 17* 2.2. This limit value corresponds to an interface tangent to the solid surface at the contact 

line above which the dynamics cannot be described by a self-similar solution anymore. It is not clear what 
the dynamics would be with a pinned contact line since direct numerical simulations of the initial set of 
equations 0 and 0 systematically lead to a finite time singularity at the boundary. We suspect that this 
singular behavior is reminiscent of the physical change observed in the experiments. When Uq is larger that 
1.5, we observe that the contact line is not pinned anymore and starts to move on the solid, as illustrated 
in the inset of figure In this case, the correct formulation would require a model for the moving contact 











6 



FIG. 4. Experimental observations of a compressed and a elongated liquid ligament, compared to the self-similar 
solution ( |16| obtained for Oh = 0.0025. The radius and capillary time of the experiment are R = 2.5 mm and 
Tc = 0.016 s, respectively. First row : The time interval between two successive frames is 0.2 ms, that is = 1/8. 
The self-similar solution (red dashed curves) has been obtained as the superimposition of two impulsive solutions, 
the impulses being located at the two red dots. The first red dot (to) corresponds to a velocity Uo = 0.645 and the 
second dot (ti) to a velocity increment Uq = 0.821. Second row : The time interval is 0.8 ms, that is St = 0.5 and 
the two dots correspond to velocity increments Uo = —0.684 and —0.489 respectively. 


line, relating the contact angle and receding speed of the contact line. 

When Uq is smaller than 1.5, the contact line remains attached and the experimental profiles can be 
compared to the self-similar solution. Such a comparison is shown in figure for a compressed and a 
elongated configuration. The experiments is done with a water ligament, hold between two vertical metal 
(brass) cylinders of diameter 5 mm. Its initial shape is slightly deformed by gravity. The upper cylinder is 
suddenly set into vertical motion and the velocity Uq is obtained by post-processing the experimental movies. 
For both cases, we find that the evolution of the velocity of the upper cylinder can be approximatively 
reproduced by superimposing two different impulses at to and ti, indicated by the red dots in figure]^ 
The theoretical solution is then obtained by adding after ti a second linear solution to the first solution 
generated at to- In both cases, this linear solution agrees quantitatively with the experimental profiles and 
the agreement is sufficiently good to confirm the relevance of the self-similar solution. These results also 
demonstrate that the short time response is barely affected by finite length effects and nonlinear effects. 
This constitutes an a-posteriori validation of the hypotheses made in the theory. 

The disappearance of the self-similar solution with a pinned contact line above a critical velocity is an 
interesting phenomenon which is worth emphasizing. In the theory, the self-similar solution disappears when 
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the tangent to the boundary becomes vertical, that is the contact angle vanishes. In experiments, we expect 
such a solution to disappear for a smaller velocity, since the contact angle cannot drop below the dynamical 
contact angle of the receding contact line. For the system water/brass/air that we have used, the dynamical 
contact angle is approximatively 30 degree, which gives using figure a critical velocity Uc ~ 1.5, a value 
consistent with the experimental observations. This, to close, underlines the importance of the boundary 
conditions in the detailed shape of the distorted ligament surface, although its space-time scaling remains 
independent of them. 
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